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Abstract. For the relativistic Boltzmann equation in Mij , this work proves the 
global existence, uniqueness, positivity, and optimal time convergence rates to 
the relativistic Maxwellian for solutions which start out sufficiently close under 
the general physical soft potential assumption proposed in 1988 I13| . 



1. Introduction and statement of the iviain results 

In the study of fast moving particles the relativistic Boltzmann equation is a fun- 
damental physical model [5l[BJ[2T] . In early work of Glassey & Strauss from 
1991, 1993, and 1995 global existence and uniqueness of nearby equilibrium solu- 
tions, and their large time convergence rates were shown in the torus (Ti^) and also 
in the whole space (Ri^). On the torus the convergence rates are exponentially fast, 
and in the whole space the convergence rates are polynomial. Their assumptions on 
the differential cross-section, a{g, 9), fell into the regime of hard potentials. Further 
results for the hard potential case can be found in [19]. However, for relativistic 
interactions, when one considers particles that are fast moving, a very important 
physical regime is the soft potential case; see [12] for a physical point of view. 
In recent work [46], the global existence, uniqueness, and rapid time convergence 
rates for nearby equilibrium solutions to the soft potential relativistic Boltzmann 
equation was shown on the torus (Ti^). However the difficult whole space case has 
remained a challenging open problem. In this work we prove the global existence, 
uniqueness, and optimal time convergence rates for nearby equilibrium solutions to 
the relativistic Boltzmann equation in M.^ under the general physical soft potential 
assumption proposed in 1988 by [l3l . 

The relativistic Boltzmann equation is given by 

(1.1) dtF+p-V,F^Q{F,F). 

The solution, F = F{t,x,p), is a function of time t G [0,oo), space cc G and 
momentum p £ K'^. It is conventional to denote the normalized velocity as 

Steady states of the relativistic Boltzmann equation are the Jiittner solutions, which 
are commonly called relativistic Maxwellians. They are given by 

(1.3) J{P) ^ 
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The collision operator, Q{F,F), is defined in (|1.4p . For the sake of simplicity but 
without loss of generality we have normalized the physical constants to be one. 
The entropy of the relativistic Boltzmann equation is physically defined as 

dx I dp F{t,x,p)\nF{t,x,p). 

Boltzmann's H- Theorem then corresponds to the formal differential inequality 

which predicts that the entropy of solutions will be non-decreasing as time passes. 
It is well known that the steady state relativistic Maxwellians (|1.3p maximize the 
entropy which grants the intuition of convergence to (|1.3p in large time. 

It is this physical reasoning that our main results, as stated below, make mathe- 
matically rigorous in the context of perturbations of the relativistic Maxwcllian for 
a general class of soft potential cross-sections in the whole space M^. 

1.1. Notation. In this section we define several notations which will be used 
throughout the article. We first introduce the center of momentum expression 
for the collision operator, as presented in |47| . In particular we have 



(1.4) Q{f,h)^ dq / d0JV,<j{g,9)[f{p')h{q')-f{p)h{q)]. 

where V0 — V0{p,q) is the M0ller velocity given by 



(1-5) W0 = U0(p,q) "^^ 



p 


q 


2 


P 




pO 


9° 




pO 



pOqO- 



Here a relativistic particle has momentum p = {p^ , p^ , p'^) G M'^, with its energy 
defined by p'^ = -^/l + |pp where |pp p ■ p. The post-colhsional momentum in 
the expression (jl.4p can then be written as: 



p = -IT- + o [^ + {i-^){p + q) 



(^g) 2 2 V''''" ' \p + q\' 

, p + q 9 ( , , -.w , dp + q)-^ 



2 2 V " ' \p + qV 

where 7 = (pP + q^)/\/s. These will satisfy (jl.7p . The angle further satisfies 
COS0 = k ■ u) with k = k{p,q) and |fc| = 1. The unit vector, fc, has a complicated 
expression as given in [47l Eq. (14)] but its precise form will be inessential. 
Now conservation of momentum and energy is given as 

(1.7) .»+,°^.;+," 

p + q = p +q ■ 

Furthermore, the relative momentum, g, is denoted 

(1.8) g ^2{pOqO-p.q-l). 
Then the quantity "s" is defined as 

(1.9) s = 2(A°-p-9 + l) > 0. 

Notice that s = g^ + 4. Also a{g,6) is the differential cross-section or scattering 
kernel; it is designed to measure the interactions between particles. We give a 
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standard warning to the reader that this notation, which is used in [6], sometimes 
differs from other authors notation by a constant factor. 

For an intcgrable function g : — > R, its Fourier transform is defined by 

r ^ 
g{k) = J'gik) / e-^''''^-''g{x)dx, x • A; = V Xjkj, k G ]R^ 

where i = \/— 1 G C is the imaginary unit. For two complex vectors a,b S C"^, 
(a I b) — a ■ b denotes the dot product over the complex field, where b is the 
ordinary complex conjugate of b. 

For a function f{t,x,p) with x £ M.^, p e R^, t G [0,oo), we define the norm 



\f{t,x,p)pdxj dpj dt 

Here G [1, oo). In the above expression, it is understood that if = c» for some 
i, then the expression is modified accordingly by replacing the L'"' norm by the L°° 
norm as usual. In particular for any G [1, oo] we have 




f\\L:^L':-'L7 ~ llll/l 



\L'-Hn)\\L'-2(R3) 



L'-i([0,oo)) 



The norm ||/||lj;2^^3 is defined similarly (and it is a function in <). The norm ||/||ij^3 
is also defined similarly (and it is analogously function in t and p). 

To study the linear and non-linear time decay rates we introduce cr^.m as 

(1.10) ^-"=1(7-^)+?' '"^o- 

Above the parameter r satisfies r e [1,2]. This <7r.m is the standard notation for 
studying time decay rates for kinetic equations in the whole space. 

We also use the norms || • Ili2(^n.) and || • |ji2(-^m) with m > 0. In this definition 

= H"^(R^) is the standard homogeneous based Sobolev space. The L^(Mp) 
inner product in the momentum variable p is denoted (•,•). Now, for £ £ R, we 
define the following weight function 

(1.11) w, = w,{p)'^{py'/'. 

The constant 6 > is defined for the soft-potentials in p.lSp below. For the soft 
potentials, we will observe later on that wi{p) ~ l/i^(p) (Lemma I2.2p . 
We also define the important temporal weight: 

(1.12) wk = wkit) - (1 + ^)^ fc > 0. 

Lastly, the notation A < B will imply that a positive constant C exists such that 
A < CB holds uniformly over the range of parameters which are present in the 
inequality and moreover that the precise magnitude of the constant is unimpor- 
tant. The notation B > A is equivalent to A < B, and A B means that both 
inequalities A < B and B < A hold simultaneously. 

Throughout this paper, furthermore C denotes some positive (generally large) 
constant and c denotes some positive (generally small) constant, where both C and 
c may take different values in different places. 
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1.2. Perturbation equation and statement of the main results. We will now 
explain in detail the main results of this paper. We define the standard perturbation 
f{t,x,p) to the relativistic Maxwellian (|1.3p as 

= J + Vlf. 

With (|1.7p we observe that the quadratic collision operator (|1.4p satisfies 

Q( J, J) = 0. 

Then the relativistic Boltzmann equation (jl.ip which will be satisfied by the per- 
turbation / — f{t,x,p) is given by 

(1.13) dtf+p-V,f + L{f)=T{f,f), f{0,x,p)^foix,p). 

The linear operator L{f) is defined in (|1.14p . And the non-linear operator r(/, /) is 
defined in (|1.17|) . They arc derived from an expansion of the relativistic Boltzmann 
collision operator (|1.4p . In particular, the linearized collision operator is given by 

(1.14) L{h) = -J-'/^Q{J, %/j/i) - J-'/^QiVlh, J) 

= i^{p)h- K{h). 

Above the multiplication operator takes the form 



(1.15) ^^(P) = / dq cL;v^a{g,9) Jiq). 
The remaining integral operator is 

Kih) / dq f doj V, a{g,e)^/jiq){^/j{^ h{p') + ^7(7) h{q')} 

(1.16) - / dq f dcuv^ a{g,e) ^J{q)J{p) h[q) 

= K2{h)~K^{h). 
The non-linear part of the collision operator is defined as 

(1.17) T{huh2) = J-^/^QiVJhi,^h2) 



dq / doJ a{g,e) ^/jiq)[hi{p')h2{q') - hi{p)h2{q)]. 

For these operators, we use the following general conditions on the kernel. 

Hypothesis on the collision kernel: For soft potentials we assume the collision 
kernel in (jl.4p satisfies the following growth/decay estimates 

<y{g,o)<g~'' MO), 

^^■^^^ a{g,e)>(^^^g-' aoiO). 

We consider angular factors < (yo{0) ^ sin' with 7 > —2. Additionally ctq^O) 
should be non-zero on a set of positive measure. We suppose <b < min(4,4-|-7). 
For hard potentials we make the assumption 

(1.19) a{g, 9) < {g^ + g-') ao{0), a{g, 0) > ( ^) g^ ao{0). 



In addition to the previous parameter ranges we consider < a < 2 + j and also 
< b < min(4, 44-7) (in this case we allow the possibility of b ~ 0). 
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This hypothesis contains the general physical assumption on the kernel which 
was introduced in |13| (and we add the corresponding necessary lower bounds). 
We are now ready to state our main result: 

Theorem 1.1. Choose £ > max{0,3/6— 1}, r E [1,6/5) and k € (l/2,(7r,o] where 
(7r,Q is given by (fTTOl) . Consider initial data fo = fo{x,p) € LlL^ n L'^{Ll fMf). 
There is an i] > such that if H/oHl^L'^ + \\we+k.fo\\L^{L^nL°°) ^ then there 
exists a unique global in time mild solution (|5.1I) . / ~ f{t,x,p), to the relativistic 
Boltzmann equation p.l3p with a soft potential kernel (jl.lSp which satisfies 

l|w'f/llL-(L2nL-)(0 < C'f,fe(l + 0"'' (II/oIIl^lj + lk'£+fe/ollL-(L2nL-)) ■ 

These solutions are continuous if it is so initially. We furthermore have the posi- 
tivity, in other words F ^ p + ,Jpf > 0, if Fq = fi + y/Jifo ^ 0. 

In the next sub-section we will discuss some historical results related to our 
main theorem to explain what has been done in the past in connection with our 
result. Explanations about why the convergence rates in Theorem 1 1 . II are said to 
be optimal can be found for instance in [101148]. 

1.3. Historical discussion. The relativistic Boltzmann equation is the primary 
model in relativistic collisional Kinetic theory. In the next few paragraphs, we will 
provide a short review of the mathematical theory of this equation. We mention a 
few books on relativistic Kinetic theory as for instance [511^1^. 

In 1988, Dudyhski and Ekiel-Jezewska [13] proved that the linear relativistic 
Boltzmann equation admits unique solutions in L^. Afterwards, Dudyhski |11] 
studied the long time and small-mean-free-path limits of these solutions. 

In the context of large data global in time weak solutions, the theory of DiPerna- 
Lions [B] renormalizcd solutions was extended to the relativistic Boltzmann equation 
in 1992 also by Dudyhski and Ekiel-Jezewska [H]. This result uses the causality of 
the relativistic Boltzmann equation [T5l[T6]. Results on the regularity of the gain 
term are given in [HIST]; the strong compactness is studied by Andreasson [1]. 
These are generalizations of Lions [37] result in the non-relativistic case. Further 
developments on renormalized weak solutions can be found in |34 | I35| . 

Notice also the studies of the Newtonian limit [HH5] for the Boltzmann equation. 
We further mention theories of unique global in time solutions with initial data that 
is near Vacuum as in Glassey [22] and [22 l l45 ] . 

Note further the study of the collision map and the pre-post collisional change 
of variables from [23]. Then [18] provides uniform L^-stability estimates for the 
relativistic Boltzmann equation. Now there is a mathematically rigorous result 
connecting the relativistic Euler equations to the relativistic Boltzmann equation 
via the Hilbert expansion as in [41] . 

We point out results on global existence of unique smooth solutions which are 
initially close to the relativistic Maxwellian for the relativistic Landau-Maxwell 
system [15], and then for the relativistic Landau [20] equation as well. Further 
[55] proves the smoothing effects for relativistic Landau-Maxwell system. And [5^ 
proves time decay rates in whole space for the relativistic Boltzmann equation (with 
certain hard potentials) and the relativistic Landau equation as well. 

Further previous results for unique strong solutions to the hard potential rela- 
tivistic Boltzmann equation are as follows. In 1993 Glassey and Strauss [24] proved 
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asymptotic stability similar to Theorem 1 1.1 1 in with £ > 3/2 in T^. They con- 
sider collisional cross-sections which satisfy (|1.19|) for the parameters b E [0, 1/2), 
a e [0, 2 - 26) and cither 7 > or 

• f„ 1 ,2-26- 
I7I < mm j 2 - a, - - 6, 



which restricts to7> —i if 6 = say. They further assume a related growth 

In |19j this growth bound was 



2 

So- 



bound on the derivative of the cross-section 

removed while the rest of the assumptions on the cross-section from |24] remained 
the same. These results also sometimes work in smoother function spaces, and 
we note that we could also include space-time regularity to our solutions spaces. 
For the Cauchy problem, under similar assumptions on the collisional cross section, 
Glassey & Strauss [25l in 1995 proved the global existence and uniqueness of nearby 
equilibrium solutions, and their large time polynomial convergence rates. 

However, it has been noted that for relativistic interactions, when one considers 
particles that are fast moving, a very important physical regime is the soft po- 
tential case [12]. Recently [46], the global existence, uniqueness, and rapid time 
convergence rates for nearby equilibrium solutions to the soft potential relativistic 
Boltzmann equation was shown on the torus (Ti^). However the difficult problem 
of determining the optimal convergence rates in the whole space case under the full 
soft potential assumption from (|1.18p remained an open problem open prior to the 
main theorem of this paper. We also believe that the methods used in this paper 
can be used to to treat the full hard potential assumption from (|1.19p , however we 
consider it to be worthwhile to carry out this extension. 

We reference other related and important previous work such as [^[51[71 [^[^[^ 
[2S{2i[30l[3Tl[33[3li0j|43|4l|4i|4^ We refer the reader to the discussions in [9] 
UnillH] for more detailed explanations of historical developments and mathematical 
methods used. We will discuss key new ideas and approaches during the course of 
the paper at the appropriate mathematical places below. 

The remainder of this article is organized as follows. In Section [5] we give the 
proof of LpL^ decay for solutions to the linearized equation (|2.5p . Then in Section 
[3] we prove the linear L^L^ time decay. Section |4] briefly explains the linear L^L"^ 
time decay. Lastly, in Section [5] we use the results from the previous sections to 
establish the non-linear time decay rates and the global existence of solutions. 



2. Linear decay theory in L^L^ 

Our main goal in this section is to prove time decay of solutions to the linearized 
relativistic Boltzmann equation with the soft potentials (|1.18p in L^L^. We consider 
the linearized equation with a microscopic source g = g{t, x,p): 

(2.1) (dJ+p.^J + Lf^g, 

I j\t=o — JO- 

For the nonlinear system (|1.13p . the non- homogeneous source takes the form of 



(2.2) 



5 = r(/,/). 
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In this case _g = {I — P}.g, is microscopic as in ()2.9p . Then it is standard to observe 
that solutions of (|2.1|) formally take the following form 

(2.3) f{t)^Uit)fo+ f dsU{t-s) g{s), t/(t) e-*(^+J^-^-). 

JO 

Here U{t) is the linear solution operator for the Cauchy problem corresponding to 
(|2.ip with 5 = 0. The main result of this section is stated as follows. 

Theorem 2.1. Fix 1 < r < 2, m > 0, and £ e M. Consider the Cauchy problem 
()2.ip with g = 0. For the soft potentials (jl.lSp with j > 0, the solution of the 
linearized homogeneous system satisfies the following time decay estimate 

(2.4) \\weU{t)foh.^H^^ < (l + t)---||^«,+5/o||L^Ls + (l + t)"^'/'lk£+,/olL.(^™), 

for any S > 2cr Recall that ar^m is given by (jl.lOp . 

We will prove this main theorem in the following few sub-sections. In Section 
12.11 we develop the theory of the different components of the solution to the linear 
equation l|2.5p . The main point of this section is to derive a collection balance law 
equations and high-order moment equations for the coefficients of the different com- 
ponents of the linear solution to the relativistic Kinetic equation. Section 12.21 then 
provides the relevant weighted instantaneous time-frequency Lyapunov inequality. 
Finally Section 12.31 contains the proof of the time decay of solutions to the linear 
equation as stated in Theorem 12.11 using the previous developments combined with 
an interpolation technique and the standard Holder and Hausdorff- Young argument. 

2.1. Linear Bounds and Decay. We consider (|2.ip with g = as 

(2.5) 9t/+p-V,/ + L(/) = 0, /(0,2;,p) = /o(a;,p), 

For the relativistic Maxwellian, J, we have the normalization J{p)dp = 1. Ini- 
tially, we introduce the notation for some integrals as follows 



/^=7 p'jdp, ^p°^'j^p^ ^^"=7 p'-^'^p^ 

JR3 



11 dcf / Pi J, ,,1122 45* / P1P2 77 ,,1111 4£f / Pi 



Jr3 (P") Jr3 (P^T 



2 

11 dcf / Pi 



Note that the constants introduced above can be expressed in terms of the Bessel 
functions but, since we have no need to explore the delicate properties of Bessel 
functions, we will only utilize the expressions above. 

We list some of the results in [46], which will be useful below. 

Lemma 2.2. 05]. Consider (jl.lSp with the soft potential kernel (ll.lSp . Then 

v{p) « (/)-''/2. 

More generally, j'^^ dq J^2 du) V0 <7{g, 9) J°'{q) ~ (p")^*"/^ for any a > 0. 
Given a small e > 0, choose a smooth cut-off function x = xid) satisfying 



(2.6) x(.9) 



1 if .9 > 2e 
if 5 < e. 



8 



R. M. STRAIN AND K. ZHU 



Now with (f2Jl) and (fTTT6| we define 

Kl'^h) f dq [ du {1~ xig)) V, a{g, 9)^7{q)^J{^) h{p') 

(2.7) +/ dq f diu (i-xig)) V0<y{g,o)^7ii)yji^hiq'). 



Define ^{h) similarly. We use the splitting K = K^->^ + K^. The following 
representation is derived in the Appendix of [15] : 



Kf{h)= / dqkf{p,q) h{q), t = 1,2. 
Here is an estimate of kf(p, q): 

Lemma 2.3. [46| . For the soft potentials (jl.lSp . the kernel enjoys the estimate 
< k^ip,q) < {p^q^y^ (p" + g")-"/' e-l''-"!, C„ c> 0, 

wii/i C =^min{2 — |7|,4 — 6, 2} /4 > 0. This estimate also holds for k^ . 

We will also use the following estimate: 
Lemma 2.4. |46j . -Fia; any small rj > 0, we may decompose K from (jl.l6p as 

K = Kc + Ks- 

Here, for any i > 0, and for some R = R{ri) > sufficiently large we have 

\{wjKchi,h2)\ < C^||1<_r/ii||l2||1<h/i2||l2, 

where 1<r is the indicator function of the ball of radius R centered at zero. We 
furthermore have the following estimate for the small part 

\{wlKshi,h2)\ < r^\\wihiv^\\^\\wih2V^\\Ll- 

For every fixed {t, x) the null space of L from (jl.l4p is given by the five dimen- 
sional space pi] : 

(2.8) = span|\/j,pi\/j,p2%/j,P3\/^,/%/j} • 

We define the orthogonal projection from L^(]Rp) onto the null space A/" by P. 
Further expand P/i as a linear combination of the basis in 



(2.9) Vh= \ a''{t,x) + Y,b){t,x)p,+c''{t,x)p''\y/j, 



where 

a" ~ I ftv jap — fi 
,h _ I^^hp^dp 

,^_ J^,hip^Vj-^i^Vj)dp 

We can then decompose f{t, x,p) as 

(2.10) / = P/ + {I-P}/. 

With this decomposition we have the coercive estimate: 
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Lemma 2.5. [46j. L >0. Lh = if and only if h = Ph. And 3do > such that 
Re(Lh,h) > (5o|k^/^{I-P}/i||?2. 

P 

We will see that for a solution to (|2.5p the coefficients of P/ satisfy balance laws: 



(2.11) dtaf (l + V. • fo/ ^) + ^ a,„A,„({I - P}/) = 0, 



3 

(2.12) at&J/i"+9,aV"+ajcV'+ I]a„,e„y({I-P}/) = 0, 1<J<3 



m — 1 

11 



(2.13) dtcf f / - j + V, • 6^ - ^ ) + E ^™A,„({I - P}/) = 0. 

\ M / \ M / 771 — 1 

These will be derived in the following developments. 

Now we are using the notation 9777 = g§~- The high order moment functions, 
Qmjih), which are included in the above balance laws, are given by 

(2.14) Q^,{h) f f ^ - a,] ,/j(p)h{p)dp, 
for 1 < m, j < 3, and ai satisfies 

(2.15) ^1122 + ^^^11 = 0. 

Our choice of ai is for a simpler expression in (|2.25p . Furthermore Am{h) is 

(2.16) A^{h) = (^-^ - ^/J(pjh{p)dp 

11 

for 1 < TO < 3, and a2 = Our choice of a2 is for a simpler expression in (|2.27p . 
We apply the decomposition (|2.9p and (|2.10p to the equation (|2.5p . Multiplying 
T5)) by PiVJ, p^\/j for 1 < i < 3 and integrating over R'^, one can get 



(2.17) dt I fVJdp+ [ p-VJ^dp = 0, 

(2.18) a*/ f^p^dp+ ( p-V.J^p^dp = Q, l<i<3, 

JR3 JR3 



(2.19) 



5t / /%/jp°dp + / P • Vccf^P°dp = 0, 

JR3 JR3 

Using (|2.10p and plugging (|2.9p into the above equation (|2.19p gives 

(2.20) dtaffi° + 9tc^/° + • 6^^" = 0. 

Plugging (|^ and (I^TTU)) into the above equation (^1^ gives (^1^ using (P?^ . 
Plugging (|2.9p into the above equation (|2.17p gives 

3 

(2.21) dta^ + dtc^pL"" + V, • bf^il^ + 9777A 777 ({I - P}/) = 0, 

711—1 

where recall A,n(/i) from (PTTB| . Then combining (P?^ and (PT^ gives (pUj) . 
Similarly, combining (|2.1ip with (|2.20p grants (|2.13p . This completes our derivation 
of the balance laws. 
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To continue, we rewrite (12.51) as 



(2.22) dtf+p- V,(P/) - -L{{I - P}/) - p ■ V,({I - P}/) = R. 

Using ([2:221) . we get 

9*8,, ({I - P}/) - e,,{R - p ■ V,(P/) - dtiPf)) 



(2.23) = e,,(i?) - 5]afc6{(/zjr _ ^^^11) _ a.6/(^mi _ c,i/,ii) 

- ata/(Aiii - ai) - 9tC-^(/i - ai/). 

From (^^2U\\ . we have 

(2.24) d,a^U^ - .,) + a.c/^^^!^^ + V. . 6/ ^"(^V"^^ = 
Combining (P?^ . ([^T^ and (j^l^ . wc get 



/ 00 

(2.25) 9te,,({i - p}/) = e,,(i?) + atc-^ ( ^(m" - "i) - + 

For j m, using (|2.22p . we get 

9te„y({i - p}/) = e„y(i? - p • v,(p/) - dt{P.f)) 



L (^-"i)^(?^.v.(p/) + at(p/)) dp 



A™(i?) 



(2.26) = e„,(i?) - (a™6j + djhDiill''^ - «i E - P>/)- 

fe=i 

Again using (j2.22p . (|2.19p and our choice of a2, we get 
5tA™({I - P}/) = A„(i? - p ■ V,(P/) - at(P/)) 

(^^ - ^2) (j5 • V,(P/) + dtiVf)) dp 

= Km{R) - dma^ifill - a2/^") - dmC^ {^1^ - a2M") 

(2.27) = A„,(i?)-9™a-^(/iiJ-a2//Ji). 
By (PH)) and (P^gj) . we have 

J2 



V 9ta™e,,({i - p}/) = f (\pI- 3ai] vldtdrni{i - P}/)dp 

= / ^VJ({I-P}9t9™/)dp. 



(2.28) 

With this calculation in mind, wc introduce the definition: 



Ir3 p" 
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On the other hand, by (|2.25|) . we have 



^at9„e,,({i-p}/) 



3 



/ , .00 

drnQniR) + id^Acf ( ^(m" - ai) - + c^i/ 



(2.29) +9„,V,.6/(a*iJ22_^iiii^ 
Combining (P?^ and (P^^ gives 



/ ,.00 

(2.30) ^a,„e,,(i?) + Zdr,AcJ ^{pl^ - ai) - /i" + ai/i" 

+ 9™V, • 6^ (/.*ii22 _ ^1111^ ^ _^ _ ^^dtd^f) . 
By p.25p and (j2.26p . for any constant /? we have 

3 

dt 5,e^.„({I _ P}/) + l3dtd^Q^^{{l - P}/) 

3 

i=i 

+ (/3 + i)d„Acf [^^^'o' - "i) - + ai/ 
+ (/?+l)9„,a„54(Mji22_^ii 



1122 



3 



(2.31) +2d„Anb{,yor - «i E E - p}/)- 

j^m k—1 



_,,1122 
POO 



Choose (3 such that 

/3 + 
3 

Then ([OT]) becomes 

3 

5t Y - P}/) + I3dtd^e^^{{l - P}/) 

3 

/,,00 N 
+ (/? + m-adtcf f ^(m" - «l) - + «l/ 

3 

(2.32) -d^drabUl]^^ ~^iYT. 9,dkAk{{I - P}/). 
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Also (|2?30l) implies 



(2.33) ^ V 9™e,,(i?) + (/3 + l)d,ndtcf i^U' - ai) - + "i/) 

^ .7 = 1 ^ 



X ' ^'^/^OO ~ 



/3 + 1 



A({i-p}ata„j). 



Combining (|2^ and (fS^ gives 



-5t 1^^^ a,e„„({i - p}/) + /39™e™„({i - p}/) 
-^5tA({i-p}a™/) 

-A ft-'^ /;"22 f) hf ,,1122 



(2.34) 



/^OO Cm Cm Om Moo 

3 3 



^ a™e,,(i?) - ^ 9je„y (i?) - /39™e™™(i?) 



j=i 

3 



i=i 



+«iEE^^^'^A'»({i-p}/)- 

j^m k—1 

We are now ready to prove the following lemma. This lemma is a key step in 
creating a time-frequency Lyapunov function later on. 

Lemma 2.6. There is a free energy functional £f red fit, k)) which is local in the 
time and frequency variables, and takes the form of 

£free{f{t,k)) ^ I ^ -^G^^ ({I _ P}/) | 



m ^ ' ' 

'15 + 1 ik, 



+«iE 



3 l + |fc|2 



^({I-P}/)I -bi 



1^ A,({I-P}7)| - 



(2.35) 



ih 



bf 



for some constant ki > 0, such that one has 

(2.36) <C||{I-P}7i.^||i. 
for any t >0 and k e M.^ . 

Proof. The proof of this Lemma 3.1 uses the overall strategy in [5], even though 
significant new difficulties arise because we are incorporating the effects of special 
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relativity. We shall make estimates on b-^ , and Hq^q-^ + /i^^c-^ individually and 
then take the proper linear combination to deduce the desired free energy inequality 
([Oe]) . Firstly, notice that 

J'a^ = a^f = a^, Fb^ = b^^ = b^, Fc^ = c^^ = J, 

and likewise for the high-order moment functions Qj„i{'), ) and A{-). 

Estimate on b^ . We claim that for any Q < 5i < /igo^^ , it holds that 



dtRcY, Y.^kjQjr,^{{l - P},f) + /3ifc™e„„„({I - P}/) 



m \ ] 



+^*fc„,A({i - p}7) I -bL ] + i^il]?'- s,)\k\^ 



bf 



(2.37) 



<<5i|fc|' 
C 



<5i|fc|^ 



-(l + |fc|^)||{I-P}/z.^||i.. 



In fact, the Fourier transform of (|2.34p gives 



^ ifcjej,n({i - p}7) + /3ifc™e™™({i - p}/) 



ik^Ai{I-P}f) 



+ i\kfbl + klbl)^^l 



122 
00 



= ai^^fc,fc,Az({I-P}/) 

+ ^ ^ ikmOjjiR) - '^ikjQjjniR) - PikmOmmiR)- 

j^m j 



Taking further the complex inner product with b(^ gives 



St ^ *%e,„,({i - p}7) + /?zfc„e„„„({i - p}7) 



f3 + l 



(2.38) 



3 

/1+/2 



*fe,„A({I-P}/) I -6" + 



6£, 



14 



R. M. STRAIN AND K. ZHU 



where 



^1 = + 

1=1 j^m 



J2 ik,e,r,r{{l - P}7) + /3«fcme™„({I - P}/) 



/3 + 1 . 



ifc„,^({I - P}/) I -dtbi 



Then Ii is bounded by 



|/i| < <^i|fc|' 



C_ 



2 C 



jm ^ — 1 

For /2, one can use the Fourier transforms of (|2.12p : 

(2.39) dthl^i^^ + ik,{J^f + cV") + - P}/) = 0' 

m 

to estimate it as 
\H<h\k\^ 



5^ 



i^iMEl©^™(ii-p}/) + ^({i-p}/) 



On the other hand, notice from (|2.22p that 

i? = -p.zHi-P}7+Mi-p}7, 

which imphes 

E|0.m(i?)f +|^(i?)f < c(i+ifcnii{i-p}Miii.. 

Thus we see that (|2.37p fohows from taking the real part of (|2.38p and plugging 
the estimates for 1\ and I2 into that. 

Estimate on . We claim that for any < 82 < Moo ~ Q^2Mo^7 it holds that 

dtReY, (Aj({I - P}./) I tk,J) + {pill - «2/iJ' - S2)\k\ 
j 

(2.40) < S2\k\' \bf\' + ^(1 + \kmi pjh^wi.. 

In fact, similarly as before, from the Fourier transform of (|2.27p 
ft A, ({I - P}7) + ik.Jifill ~ a2fil') = UR), 
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one can get 

dt (a,({I - P}/) I ik,J) + |fc,f \af^ {fill - «2m") 
(2.41) = (a,{R) I ik,J) + (a, ({I - P}7) I ik.dtJ) =h + h. 



Now is bounded by 



l^3| < 52\k, 



2 C 



and from the Fourier transform of p. lip 

_,,11,,00 _ ,,11 ,,0 ^ 

d^-' + • ^'^^Sr^^ + E ^fc. A, ({I - P}/) 



= 0, 



where /4 is bounded by 

|/4| < 52\k\'^ 



hi 



Notice that similar to ,„i, it holds that 



i||2 



A,(i?) <c(i + ifcn!i{i-p}/z.^ii 



LI- 



Then, ()2.40[) follows from p.4ip by taking summation over i, taking the real part 
and then applying the estimates of /a and /4 . 

Estimate on fJ-l^a-^ + fi^^c-^ . We notice that c-^ is a linear combination of iJ-l^a^ + 
fj}^c^ and a-^. So our estimates on nl^a-^ + fi^^c-^ and a-^ imply the estimate of c-^ . 
We claim that for any < < it holds that 



(2.42) a^RcE (bl I ikmifil^af + fi^'J) 



63 \W 



6/ 



'+^ifcni{i-p}Miii,. 
03 p 



In fact, by taking the complex inner product with ikj{fil^af + fi^^cJ) and then 
taking summation over 1 < j < it follows from (|2.39p that 

5,E(^^Mfc.(/^JV + /V)) + '^''' 

j 



/II 



■e,™({i-p}/) I ik.ifil'J+fi'^J) 



= E 

J 

(2.43) =/5 + /6. 

Now /s is bounded similar to previous estimates as 

\h\<S3\kf 



' + ^ifcni{i-p}Miii.. 

03 p 
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For Iq, it holds that 

\i6\<c\k\' b^\c\km{i-p}M\\i,, 

where we have used the Fourier transform of (|2.1ip as 

_,,ii,,oo _ ,,ii,,o ^ ,,00 

and the Fourier transform of (|2.13p 

Therefore, (|2.36p fohows from the proper linear combination of (|2.37p . (|2.40p and 
(|2.42p by taking < 61,52,63 < 1 small enough and also kq > large enough. This 
completes the proof of Lemma 12.61 □ 



2.2. Weighted time-frequency Lyapunov inequality. In this subsection, we 
shall construct the desired time-frequency Lyapunov functional. 

2.2.1. Estimate on the microscopic dissipation. The first step in our construction of 
the time-frequency Lyapunov functional is to estimate the microscopic dissipation 
on the basis of the coercivity property in Lemma 12.51 of L. 

Consider (|2.5p . taking the Fourier transform in x grants us 

(2.44) dtf + ip- k f + Lf ^0. 

Then we multiply equation (|2.44[) with f{t, k) and integrate over to achieve 

iaj/||2,+Rc(L/,/)-0. 
By Lemma [^751 one has that 

(2.45) dml.+\\\u^'^{l-V]f\\^<Q. 

This is the first main estimate which we will use in the following. 

2.2.2. Macroscopic time- frequency weighted inequality. In this section we prove the 
following instantaneous Lyapunov inequality with a velocity weight £ G R: 

(2.46) P}./lli. + \h"^w,{l - P}/Ili. 

<Cxmu^'^f\\l.+C\\l<c{l~'P}f\\%. 

We split the solution / to equation (|2.5p into / — P/ + {I — P}/, take the Fourier 
transform as in (|2.44p . and then apply {I — P} to the resulting equation: 

dt{l - P}/ + ip • k{l - P}/ + L{1 - P}/ 

= -{I - P}(ip • kVf) + P(ip • fc{I - P}/). 

Multiply the last equation by W2i{'i. — P}/ and integrate in to obtain 

(2.47) ~\\w,{l-V}f\\l,+Rc{w2iL{l-V}f,{l-V}f)=T^, 
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where 

Ti = -Re ({I - P}{ip ■ fcP/), W2i{I - P}/ 

+ Re (P{ip ■ k{I - P}/), W2di - P}/ 
As a result of the rapid decay in the coefficients of p.9p we obtain 



iril < v\W'^'-wdi^P}f\\h + c,\kf (||u,_,{i-P}/||2 + IW'^'PfWh) , 



which holds for any small ?/ > and any large j > 0. 
By (|1.14p and Lemma [2131 we have 



Re{w2eL{I - P}/, {I - P}/) > A||u.,{I - P}.^ 



4l|2 



q|i<c{i-P}/||i.. 



This holds for a small A > and a large C > 0. (We have used the fact that Lemma 
12.41 indeed holds for any ^ e R, as can be seen from the proofs in [JHl Lemma 3.3] 
and [46l Lemma 3.6].) Plugging the last few estimates into (|2.47p proves (|2.46p . 
We furthermore remark, following the same procedure as above, that we obtain 

(2.48) + MW'^'^JWli < CWUcfWl.- 

In other words, if we multiply (|2.44p by w'^^f{t,k), integrate in M.^ and use the 
same estimates as in the last case it follows that we obtain (j2.48p . 

2.2.3. Derivation of time-frequency Lyapunov inequality. Now we prove 

Theorem 2.7. Fix £ e M. Let f be the solution to the Cauchy problem (|2.H) with 
g — 0. Then there is a time-frequency functional £({t, k) such that 



(2.49) 



£iit,k) 



where for any t > and fc G M'^ we have 

(2.50) dMt, fc) + A (1 A |fc|2) Wu^'^wJWl. < 0. 

77ere 1 A |fc|2 'J min{l, |fc|2}. 
Proof. We first define 

(2.51) £{t,k)''^'\\f\\l.+K3£free{t,k), 

for a constant K3 > to be determined later, where £free{tik) is given in (|2.35p . 

One can fix K3 > small enough such that £{t, k) « |!/||?2- 

_____ p 

A linear combination of (|2.45p and (|2.36p implies that 

(2.52) dMt, k) + A||j.i/2{I - V}f\\l, + 



+ \c?] <0, 



where note further that one has \a\ 



< 



|p/lli.- 



To do the weighted estimates, in particular for the soft potentials (jl.lSp . we need 
to use the energy splitting as follows. With ()2.5ip we define 



£^it,k) 
£}{t,k) 



L|fc|<l 



L|fc|>l 



(^£it,k) + K4we{I-F}f\\ 



£{t,k) + K5\\wif\ 



18 



R. M. STRAIN AND K. ZHU 



Again K4, K5 > will be determined later. 

We prove estimates for each of these individually. For £({t, k) we combine ()2.52|) 
with (|2.48p for |fc| > 1 to obtain for a suitably small > that 

Here we have used the fact that when |fc| > 1 then j^j^'^^^p > i. 

Furthermore, when |fc| < 1 it holds that j^j^j"^^^,|2 > In this case we combine 
([2?52| with (|2^ on |fc| < 1 to obtain for a small K4 > that 

dt£^{t,k) + A|fc|2|ii.i/V/||ia|fe|<i < 0. 

Lastly we define £e{t, k) Ei{t, k) + El{t, k) and we notice that (|2.49p is satisfied. 
Then (|2.50p follows from adding the previous two differential inequalities. □ 

2.3. Proof of time-decay of linear solutions. Our proof of Theorem 12.11 is 
based on Theorem 1 2 . 71 and the interpolation argument given below. 

Proof of TheoremliJi We define p{k) = A (l A \k\^) . By ({2301) we have that 

£e{t,k) <£e{0,k), 

for any £ G M. Now we use the interpolation technique as in [33], but in a different 
context. In particular, for j > 0, using (|2.49p we have 

£dt,k)<£iL^^^'\t,k) £ji^^'^(t,k) < 

We therefore conclude that 

£l'^'^^'{t,k) < \W'/'wJ\\l.£ji]it,k) < \W'/^wef\\h/,i]{0,k). 

Now we can rewrite (|2.50p . for any A: S R'^, as 

dMt, k) + xp{k)^i^^^'\t, k)£;^j\o, k) < 0. 

To prove (|2.4p . one can bound £{t, k) as follows 

dMt,k)£-'-^''{t,k) < -p{k)£;^^'{o,k). 

Integrating this over time, we obtain 

j£;'/'{0, k) - j£i''\t, k) < -tpik)S7+fiO, k). 
For any ^ G M and j > 0, uniformly in A: e M"^. we have shown that 

£eit,k)<£e+,{0,k) + l) 

We also just used the estimate £i{0, k) < £ij^j{Q, k). 

As before, we integrate over k and split into |fc| < 1 and |fc| > 1 to achieve 

/ dk \k\'^'^£i{t,k) <(i + l\ ^ [ dk \k\^"'£^+j{Q,k). 

J\k\>l \J J J\k\>l 
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Alternatively, when \k\ < 1 we choose j to be any number j = S > 2ar.m and obtain 

dk \k\^"'£e+siO,k) ( 

'|fe|<i ■J\k\<i 



dk \k\^"'£e{t, k)< I dk |fc|2"£f+5(0, k) f ^ + 

J\k\<i \ J 



<(i + i)-^'^-||^,+,/olli.L,. 

For 1 < r < 2, this uses the same Holder and Hausdorff- Young argument. □ 

3. Linear decay theory in L'^LI 

Now we work on the linear L^L^ bounds and time decay. As is customary, we 
express solutions, f{t,x,p), to (|2.5p with the semigroup U{t) as 

(3.1) f{t,x,p)^{U{t)fo}ix,p), 
with initial data given by 

{U{0)fo}{x,p)^fo{x,p). 
Then we have the following result. 

Theorem 3.1. Fix £ > 0, r e [1,2] and k e [0, (Tr,o]- Suppose we+kfo e L^Ll, 
then under (|1.18[) the semi-group satisfies the estimate 

\\wdU{t)fo}\\L^Li < Cil + t)-" [Wtve+kfoh^Li + II/o||l|ls) . 

Above the positive constant C — Cg,k only depends on I and k. 

A key idea in the proof is that, instead of placing everything in the L^L^ space 
as in the proof of Theorem 4.1 in [46], we use the L^L^ space. The crucial new 
element that we now use is Minkowski's inequality in the following form 

\\fiy)\\Lz'^('^y)^ i<_p<oo, 

where Q is any measure space with sigma-finite measure v. This allows for substan- 
tial simplifications over previous work. Also since the inequalities used to deal with 
the term ^^""''^ (which will be defined in the following discussion) in the proof of 
[46l Theorem 4.1] do not work well in our setting, and we use a completely different 
approach. The term ^^""''^ -^^3,3 estimated with a complicated change of variables 
in the proof of [46[ Theorem 4.1]. Our new approach does not need to use any such 
change of variables as a result of our utilization of Minkowski's inequality. 

Now we first consider solutions to the linearization of (|2.5p with the compact 
operator K removed from (|2.5|) . This equation is given by 

(3.2) (at+p-V, + i/(p))/ = 0, fiO,x,p)^fo{x,p). 

Let the semigroup G(t)/o denote the solution to this system p.2|) . Explicitly 
G{t)fo{x,p)''^'e-''^P^'fo{x~pt,p). 



[ fivHdy) 




Jn 
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For the solution {U (t) fo} {x , p) , by iterating twice (as did Vidav [50]) we have 



{U{t)fo} ix,p) = G{t)fo{x,p) + I dsi Git - si)K'-^ {Uisi)h} ix,p) 







(3.3) +/ dsi G{t~si)K^G{si)fo{x,p) 





dsi f ' ds2 G{t - si)K^G{si - .S2)A-i-x {U{s2)fo} {x,p) 
Jo 

dsi / ' ds2 G{t - si)K^G{si - S2)K^ {Uis2)fo} (^,p)- 

JO 

Equivalently 

{U{t)fo}{x,p) = Hi{t,x,p) + H2{t,x,p) + H3{t,x,p) + H4{t,x,p) + H5{,t,x,p), 
where 

H,it,x,p) er^^^^'Ux~pt,p), 

H2it,x,p) f ds,e-'^P^^'~'^^K'-^{U{s^)fo}{yi,p), 

Jo 

H3{t,x,p) '=^' / dsi e-''(f)(*-^i) / dqi fc>^(p,gi) e-''(«i'-'Vo(yi-<7isi,gi). 
Just above and below we will be using the following short hand notation 

yi "= x-p{t-si), 
(3.4) 2/2 yi - qi{si - S2) ^ X ~ p{t ~ si) ~ qi{si - S2). 



We are also using the notation = \/l + and 91 = {qii , qi2 , qis) & with 
91 ~ li/li- Furthermore the next term is 

H^{t,x,p)= f dqi k^{p,qi) [ dsi r ds2 e-''(P)(*-^i)e-''('?i)(«i-«2) 
J«? Jo Jo 

xifi-x{C/(s2)/o}(y2,<Zi). 
Lastly, we may also expand out the fifth component as 

(3.5) H5{t,x,p)^ f dqik^ip.qi) f dq2k^{qi,q2) f dsi e-^^^P'^'-^^) 
Jr3 jTsia Jo 

X r ds2 e-''('i)(-^^-^^){[/(s2)/o}(y2,<Z2). 



Jo 

We have the following estimates for the Hi terms above (1 < i < 5). 
Lemma 3.2. Given £ > 0, for any k > we have 

\\wiHi\\L^L2 + \\weH3\\L^L2 < Ci^kil +ty''\\wk+efo\\L^Ll- 
Lemma 3.3. Fix £ > 0. For any small r] > and any k > we have 
\\wiH2\\L^Ll + \\weHi\\L^Ll < v{l + ty'^WwkWefWL^^Ll- 
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Lemma 3.4. Fix t > Q, choose any (possibly large) j > 0. For any small rj > 
and any k > we have the estimate 

Jo 

By the L^L^ decay theory from Theorem \2.1\ and Provosition \3.5\ (below), we can 
choose j > niax{2cr,.,07 6/6 + k} and take k e [0, cr,-,o] to have 

j\s e-^(''^^\\w^,fU.L.{s) 

< + t)-''\\w2k~j.fo\\LlLl + C„il + i)-'!|/o|U^LJ 

< + t)-' {\\we+kM\L'^Ll + ll/olUjLj) . 

The above estimates hold for any k € [0, (Tr,o] md £ > 0. On the other hand, for 
the last term involving Ri{f) if we restrict k £ [0, 1] then V?7 > we have 

\\mRiif)\\L^Ll < vi'^ + t)-''\\wkWef\\L^^L2. 

This term Ri is defined in liS.lS]} during the course of proof. 

Combining Lemma [3?2l Lemma [331 and Lemma [33] gives Theorem l3.1l Actually, 
Lemma 13.21 Lemma 13.31 and Lemma 13.41 imply that for any 77 > and k € [0, (t,.,o] 
we have 

Wmfh^Llit) < Cf^^kA'^ + ty'^Wive+kfoh^Ll 

+ 77(1 + O^'ll W/||lj^,L^ + + t)-'||/o||LJLj . 

Equivalently 

\\''^kW^f\\L^^Ll < Ce^k.jjWwe+kfoWL^Ll + C'r;||/o||L2L!;;- 

With this inequality, we have proved Theorem 13 . 1 1 sub i ect to Lemmas 13.21 13.31 and 
13.41 We now prove those lemmas. 

We will use the following known basic decay estimate, as in [46l Proposition 4.5]: 
Proposition 3.5. Suppose without loss of generality that A > /i > 0. Then 

Jo {l + t-s)>^{l + sy - {l + t)P' 
where p — p{\, fi) — min{A + /i — 1, /i} and 

0<CA,^(t)-c| iog(2 + <) 'Ifltl. 
Furthermore, we will use the following basic estimate from the Calculus 
(3.6) e"''^(l + y)'= < max{l,e''-'^fc'^a-'^'}, a,y,k>0. 

We are now ready to prove the lemmas above. 
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Proof of Lemma \3/A We first look at Hi . The following result is shown in the 
proof of Lemma 4.2 in [46]: 

(3.7) e-'^(f)* < Ckpf'\l + t)-'^ < Cuwk{p){l + Vi, fe > 0. 

So we have 

\wt{p)Hi{t,x,p)\ = \wt{p) e~''<^P^'foix-pt,p) 

< Ck \we{p)wk{p){l + t)-'' foix ~ pt,p)\ . 

So 

\\wiHi\\l2 < Ck\\wiWkZu^kfo\\Ll- 

Thus 

IhML^Li < Ce,k{l + ty'Wwk+ifoh^Li- 
Now wc turn to H3. Since 

g-i/(p)(t-si)g-iy(g)si ^ g-iy(max{|p|,|g|})t^ 

where :^(max{|p|, \q\}) is evaluated at max{|p|, \q\}, we have 
\we{p)H3{t,x,p)\ 

<Wi{p) f dsie-^^-P^' [ dqi\k^ip,qi)\\foix-p{t~Si)-qiSi,qi)\ 
Jo "'|p|>ki| 

+ w,{p) f dsi f dqi |fc>^(p,gi)| e-'^(«i)*|/o(a;-p(t-si)-giSi,gi)|. 
Jo "'ipKkil 

So by Minkowski's inequality, we have 

WwiHsWLi < Wi{p) f dsi e-^'ff)* / dqi \k^{p,qi)\ \\fo\\Ll{qi) 

Jo "'|p|>l9l| 

+ wt{p) I dsi [ dqi \k^{p,qi)\ e-^^^^^'llfohliqi). 
Jo "'ipKkil 

We will estimate the second term, and we remark that the first term can be handled 
in exactly the same way. As in Lemma 4.2 in [46j, and similar to (j3.7p . we have 

e-'''^^'*ll/o||L?(gi) <Cfc(i + t)-'-'w^fc+i(gi)||/o||L?(9i). 



Next we use the estimate for {p, q) from Lemma 12.31 When using this estimate 
we may suppose \qi\ < 2\p\. For otherwise, if say \qi\ > 2|p|, then we have |p— gi| > 
|gi|/2 which leads directly to 

(3.8) u;fe+i(gi)u;,(p)e-^lf-9il/2 < Cwk+i+i{qi)e-'^'''\/^ < C. 

In this case we easily obtain the following estimate: 

dsi Wiip) f dqi 1|,,|>2|,| \k^ip,qi)\ e-^^^^^'Wfohliqi) 

J\p\<\qi\ 

< (l+t)-^-||«;_,/oi|L=fL?, Vfc > 0, J > 0. 
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Thus in the foUowmg we assume \p\ < \qi\ < 2\p\. On this region we have 



(3.9) 



dsi we{p) 
t 



bl<kil<2|p| 



k+l 



'\p\<\qi\<2\p\ 

< Ck{l +t)-''\\we+kfo\\L^Ll "'i(p) 

'|p|<|gi|<2|p| 

Then from Lemma 12.31 we clearly have the following bound 



dqi we{qi) Wk+iiqi) \k^{p,qi)\ H/ollL^lgi) 

dqi \k^{p,qi)\. 



wi{p) 



dqi \k^{p,qi)\<C. 



'|p|<l9i|<2|p| 

This completes the time decay estimate for H3, and our proof of the lemma. 



□ 



Proof of Lemma \3.3\ As in the proof of Lemma 4.3 in [55], we will use the following 
lemma (Lemma 4.6 in |46|). 

Lemma 3.6 (3^). Fix any € > and any j > 0. Then given any small i] > 0, 
which depends upon x in (|2.6I) . the following estimate holds 

w,{p)Kl-''{h)ip)\ < rje-^P°\\w^,hU^. 

Above the constant c > is independent of rj and i — 1,2. 

From (|2.7p and Minkowski's inequality, we have 

\\Kl-^ih)ip)U2<Kl-'^{\\h\\Ll)ip)- 
So with the lemma above, we have 

(3.10) \\w,{p)Kl-^{h){p)\\r^. < ^e-^'^"\\w^,hU^r.l- 

Notice that in the result above, we use ||u;_j/i||ioo^2 instead of ||w_j/i|jioo^2 to 
avoid confusion, since p is a variable in the left- hand- side. From now on, we always 
use q to denote the momentum variable in the norm expression whenever p has 
been used to denote a particular momentum value in an inequality. For H2, with 
the result above, for any small 77' > we have 



\weH2\\Llit,p) = Wi{p) 



ds, e-'^(rf(*-^i)A'i->^ {U{si)fo} {yi,p) 



LI 



<r,' e-^P \M\\L^,Li 



dsi e-'^(P'(*-"l)(l-KSl)-^ Vj>0. 



For any A > max{l, k} we have 

\\weH2\\Li{t,p) 

<r,' wx{p)e-'P°\\mkf\\L^^Ll [ ds {I + t - s)-\l + s)-'' 



< 77 (1-^0^1 Nfe/llL-,L^ 
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which fohows from Proposition [3?5l This is the desired estimate for H2. For H4 we 
once again use p.lOp . for any small 77' > 0, to obtain 



\\weH4L2{t,p) < [ dq, k>^{p,q,) f ds, r ds2 e-^^f^^'-^i) 

Jo Jo 

<v'\M\\l^,li [ dq^ k^{p,q^)^^e-^'^° 

X f dsi r rfs2 e-''(P)(*~"i)e-''('?i)("i-"^)(l + S2)-^ 
Jo Jo 

For the time decay, from Proposition 13. 5[ we notice that 
Jo Jo 

< wx{p)wx{qi) f dsi {1 + t- si)-^ f ds2 (1 + Si - S2)-\l + 82)-^ 
Jo Jo 

< Wx{p)wx{qi){l + t)-''. 

Above we have taken A > max{l, k}. Combining these estimates yields 
\\Mp)H4\\lI 

<ri'{l + t)-''\\-cukW(f\\L^^L2 / dqi k^ip,qi)wi+x{p)w^i+x{qi)e-'='i". 



To estimate the remaining integral and weights we split into three cases. If either 
2I91I < \p\, or \qi\ > 2\p\, then we bound all the weights and the remaining momen- 
tum integral by a constant as in p.8p . Alternatively if ^\qi \ < \p\ < 2|gi|, then the 
desired estimate is obvious since we have strong exponential decay in both p and 
qi . In either of these cases we have the estimate for H4 . □ 



Proof of Lemma \3.4\ We will utilize rather extensively the estimate for k^ from 
Lemma [231 We now further split H^{t,x,p) as 

(3.11) H^{t,x,p) ^ H'^'''\t,x,p) + H'r{t,x,p), 

and estimate each term on the right individually. For M » 1 we define 

(3.12) l/iig/i = l|p|>A/l|gi|<Jl/ + l|gi|>A/- 

Notice l/iig/i + l|p|<j\/l|(ji|<j\/ = 1. Now the first term in the expansion is 

H^'<^\t,X,p) = f dqi k^{p,qi) f dq2 k^{qi,q2) lh^gh f dsi e-^^^^^'-'^^ 

JR3 Jo 

ds2 e-^(«^)(^^-^^UC^(s2)/o}(y2,<?2). 
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The proof of Lemma 4.4 in [4Q, and also p.7p . shows that for any A > we have 

ft ysi 



Jo Jo 



ds2 e 



-l^(qi)(si-S2) 



<Cxwx{p)wx{qi) dsi / ds2{l + {t-si))-^{l + {si-S2)y 
Jo Jo 

When cither |p| > M or \qi \ > M, by Lemma [^751 we have the bound 

qi)\< CM-^ + gj)"'/' e--IP-9i|. 
If either \p\ > 2\qi \ or \qi\ > 2\p\ then as in p.Sp we have 

Thus by combining the last few estimates we have 

wi{p) / dqi |fcX(p,(7i)| / dq2 \k'^{qi,q2)\ Ihigh (1|p|>2|9i| + l|p|<i|qi | ) 

Jo Jo 



< 



^||z..z../|U,,^., / d,ie---- / d,2e 



c|p-9l| / An„ p-c|gi-g2| 



;* , /^^ , (1 + S2)-'' 

X / dsi I ds2 



Jo 



(l + (t-Sl))^(l + (si-.S2))^- 



With Proposition 13.51 for any fc > and A > max{/c, 1} the previous term is 
bounded from above by 

This is the desired estimate for M ^ 1 chosen sufficiently large. We now consider 
the remaining part of H^''^^ . As in the previous estimates and p.Sp . if either 
I92I > 2|gi| or \qi\ > 2\q2\ then for any > we have 



we{p) dqi \k^{p,qi)\ / dq2 \k^{qi,q2)\ (lk2|>2|gi| + 1|,,|< 1 ) 
JK3 -/MS ^ ^ 

X HlPl<k.|<2|p|l«.'> /* ds, ds2 e-^('?^)(^-^^) 



10 JQ 

X \\{U{s2)fo}{q2)\\Ll 



Jo 



X / dsi I ds2 ~^ 



(l + (f-Sl))^(l+(si-S2))^ 
< 



Above we have used exactly the same estimates as in the prior case. Both of the 
last two terms have a suitably small constant in front if M is sufficiently large. 
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j-high 



Thus the remaining part of i?5*^ to estimate is Ri{f){t) which is defined by 

Rl{f){t)= dqi k^{p,qi) / ^92 ^^(91,92) ll|p|<|gi|<2|p|li|gi|<|g2|<2|gi| 

(3.f3) xU,gh f dsi e-'-^fH^-^i) r e-'^(«^)(^i-^=){C/(s2)/o}(y2,'?2). 



Since all the momentum variables are comparable, we have 
\\Rl{fm\\Ll< f dq, \k^p,q^)\ [ dq2 \k^{quq2)\ 



^ l|bl<l9l|<2b|l||9i|<|g2|<2|gi|l'i«Sft. 



X / dsi e-^''(?i)(*-^i) / ds2 e"^''(9l)(^■l-^^)||{C/(s2)/o}(y2,g2)|lL^ 



Next using similar techniques as in the previous two estimates, including (j3.7p . we 
obtain the following upper bound for any k e [0, 1]: 

we{p)\\R,{m)h. < ^\\r^,wd\\L-Ll f dq, (90)-''/2-Ce-|p-,,| 



X / dq^ {qlybn-C^-c\q^-q,\ ^2^25(<Zl) 

ds2 



dsi 



{l + {t-s,)Y+^ 



(1 + (S1-S2))l+^(1 + S2)fc 

< 



In this computation we have chosen (5 > to satisfy <^ < C where C > is defined in 
the statement of Lemma l2.3l This choice guarantees that W2+2i5('i'i)(<Zi)"''^^^ < C. 
We are ready to define the second term in our splitting of . It must be 

^^5°'"(*>a;,p) = l|p|<M / dqik^ip,qi) / dq2 k^{qi,q2) 

J\qi\<M "'k^ 

X /* ds, e-'-'^P^^'-'^^ r ds2 e-'^(''^)(^^-^^){C/(s2)/o}(y2,(?2). 



For any small k > 0, we further split this term into two terms, one of which is 

Hi°^^^it,x,p)^^' [ dq,k^{p,q,) [ dq2k^iqi,q2) P ds^ e^^'^P^^*-^^^ 
J\qi\<M Jr^ Jo 

xl|p|<M r ds2 e-'^('?iH^i-^=){[/(s2)/o}(y2,92) 
Jo 

+ [ dqik^{p,qi) f dq2k^{q^,q2) f ds, e-'^P^^'-''^ 

J\qi\<M JR3 Jk 

X l|p|<M r '^'^ e-'^'^^^^^^~^^-HUis2)fo}iy2,q2). 



Si — K 



The other term in this latest splitting is defined just below as iJg"™'^. Since p and 
qi are both bounded by M, from Lemma 12.21 we have 

(3.14) 1|p|<m1|,i|<a/ e-KP)(*-«i)-^'to)(^i-^2) < ^-Cit-s.)/M^^\ 
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Then for the first term in ^^"■"•'^ above multiphed by 'we{p) we have the bound 



w,{p) dq, \kHp,q,)\ / dg2 1^:^91,92)1 / d5i e-''(f)(*-^i) 

J\qi\<M Jr'^ Jo 







Si 



Jo Jo 

We obtain the desired estimate for the above terms by first choosing M large, and 
second choosing k = k{M) > sufficiently small. 

For the second term in multiplied by wiijp) for any fc > we have 



wi{p) I dq, \k^ip,qi)\ f dq2 |fc>^(<Zi, 92)! f ds^ er^^P^^'-'') 

J\qi\<M Js? Jk 

X l|p|<M ds2 e-'^(«^)(^^-^=)|i {C/(s2)/o} iq2)\\Ll 

J Si —K 



< CM\\^kWif\\L-^^Ll 

X f ds, r ds2 e-C(*-^i)/^^''%-C^(-^i-«^)/^^'''(f + 52)-^ 



Si—K 



Since S2 G [si — k,si] and k G (0, 1/2), wc know (1 + S2) > (| + si) • Then the 
previous display is further bounded above as 



< CMK\\vukWef\\L^^Ll / dsi e-^(*-^i)/^^ (1 + s.r'^ 

< Cmk(1 + ty''\\wkWif\\L^^L2. 

In the last step we have used Proposition 13.51 We conclude the desired estimate 
for by first choosing M large, and then k > sufficiently small. 

The only remaining part of Hl°^{t^ x,p) to be estimated is given by 

H'°^'^{t,x,p)''^' [ dq,k^{p,q,) f dq2k^{q,,q2) f ds, e-''^P'>^'-'^'> 
J\qi\<M JR3 Jk 

xlbl<M r " ^^^2 e-'^('i)(-^^-^^){[/(s2)/o}(y2,<?2), 
Jo 
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Now for H^^°^'^{t,x,p), we notice using Cauchy-Schwartz tliat 



ki|<M 



dqik^{p,qi) / dq2 k^{qi,q2) \\{U{s2)fo}\\Ll(q2) 



< I dqi k^(p,qi) 

'\qi\<M 



dq2 \wj{q2)k^{qi,q2)\' 



1/2 



X / dq2 



W.,{q2)\\{U{s2)fo}\\LlAq2) 



1/2N 



< CM\\W-jf\\L2 js2). 

Then 



These estimates hold Vj > 0. Above we used the sphtting as in 
furthermore, as in the previous estimates inchiding (|3.14l) . we have 



\\weHi"-^'h^Ll < Cm f ds, r ^ ds2 e-C(*-«i)/M^'%-c(.i-s.)/A/V= 

X \\u'-j{U{s2)fo}\\Ll^ 

< Cm f ds, e-*(*--)/^^'^^ /* ds2 e-#(*--)A^^^%-* (---)/^^'^^ 
Jo Jo 

X ||u;_j{C/(s2)/o}||l2^. 

Notice that the first exponential controls the si time integral, and the second and 
third exponential control the remaining time integral as follows 

<Cm f ds2 e-*(*-^^)/*^'''|k-,{C/(s2)/o}|U. . 



□ 



This completes the proof of Lemma [3^ by first choosing M large. 



4. Linear decay theory in L'^L'^ 

Now we work on the linear L^L^ bounds and decay. The argument that we 
will use in this short section is not new. By Lemma 4.2, Lemma 4.3 and Lemma 
4.4 from [46], for ^ > 0, fc > 0, small 77 and (possibly large) j > 0, we have 

\\wdU{t)fo}\\L^L^ 

< Cil + t)-''\\we+kfo\\L^Lr + ^(1 + ty'^W^kWi/WL^^LS' 

+ C, fds e-"(*-^)||7«_,/|U. (s). 



Note that even though the lemmas just cited from were proven in the context 
of Tj^ , they generalize directly to Ri^ without modification. 

Using the same upper bound for ds e~''*^*~*' || w_j/||i2 (s) as in Lemma [3.41 
we then have the following result. 
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Theorem 4.1. Given ^ > 0, r G [1,2] and k G [0, (Tr,o]- Suppose that initially we 
have wg-^^kfo G (L^ H -/j!^) , then under (|1.18[) the semi-group satisfies 

\\wi{U{t)fo}\\L^Lr < Cil + ty^ {\\wi+kfo\\L^(LlnL^) + 1I/o1Iljls) . 
Above the positive constant C = Ce,k only depends on £ and k. 

5. Nonlinear decay theory and global existence 

Suppose / = f{t,x,p) solves (|1.13p with initial eondition f(0,x^p) ~ fo{x,p). 
We may express mild solutions to this problem p.lSp in the form 

(5.1) fit,x,p) = {Uit)fo}ix,p)+N[fJ]it,x,p), 

where we have used the notation 

N[h,hKt,x,p)''^' f ds{U{t-s)T[h{s),h{smx,p). 
Jo 

Here as usual U{t) is the semi-group p.l[) which represents mild solutions to the 
linear problem (|2.5p . The purpose of this section is to prove Theorem ll.il 
As a first step, we will use the following non-linear estimate. 

Lemma 5.1. Considering the non-linear operator defined in (jl.l7p with (jl.lSp . we 
have the following pointwise estimates 

\\weT{hi, h2)\\LinL^{p) ^ '^ip)\\wehi\\L^(^LlnL^)\\weh2\\L-^(LlnL^)- 
These hold for any i > 0. Furthermore, 

||w£+ir(/li,/l2)||too(iinL~) < |lw£/li|lioo(i2nL~)||w£/l2|lL~(L2nL~)- 

Note that just as in the previous section, we use q to denote the momentum 
variable in \\wihi\\L,^(^i^2(-^i^^-)\\wih2\\L'^(LlnL^) to avoid possible confusion, since p 
is a variable in this inequality. The lemma above combined with Proposition 13.51 
will be important tools in our proof of Theorem 11.11 

Proof of Lemma[5Ji We recall pTTS]) . ([TTT)) . and pTTTj) . For £ > 0, it follows from 
(HHl) that 

m{p) < {pT" < {p'y"\qy"^ < we{p')w,{q'). 
A proof of this estimate above was given in [211 Lemma 2.2]. Thus 

W£(p)||r(/li,/l2)|iLinL- 

< 



dujdq Vf) <j{g,0) y^J{q) we{p')we{q')\\hi\\L2 L^{p')\\h2\\LlnL^{q') 



dx^dq V0 (7{g,e) V J il) W£{p)\\hi\\LinL^ (p) II 1| L^nL^ (?) 

< \\wihi\\L^^^L2r,L^)\\wth2\\L^(^LlnL^) / dujdq (T{g,0) J^''^(q) 

< zy(p)||w^/ii||i=o(i2nLoo)||w^/i2||i=o(i2nL~). 

The last inequality above follows directly from Lemma 12.21 since both the inte- 
gral and vip) have the same asymptotic behavior at infinity. That yields the first 
estimate. For the second estimate we notice from the first estimate that 

\\wi+iT{hi,h2)\\LinL^{p) ^ Wi{p)v{p)\\wthi\\L^^L2^^L^)\\wih2\\L'^(^LlnL^)- 
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But wi{p)v{jj) < 1 from Lemma [2.21 and This completes the proof. □ 

Proof of Theorem Wc wiU prove Theorem 11.11 in three steps. The first step 
gives existence, uniqueness and time decay via the contraction mapping argument. 
The second step will establish continuity, and the last step shows positivity. 

Step 1. Existence and Uniqueness. When proving existence of mild solu- 
tions to ()5.ip it is natural to consider the mapping 



M[f] {U{t)f„}{x,p) + N[f, f]{t, x,p). 
We will show that this is a contraction mapping on the space 

M^:^ = {/ : [0, oo) X X ^ R I Ww^wefh^^^LlnL^) < R}, R> 0. 

Here we recall the temporal weight p.l2p . Wc first estimate the non- linear term 
N[f,f] defined in the equation display below (jS.ip . We apply Theorem 13.11 and 
Theorem 14. 11 with £>0, and k G (l/2,(Tr.o] to obtain 

we{p)\\N[f,j2]\\LinLrit^P) 



< 



ds w,{p)\\{U{t-s)T[f,{s)j2{s)]}hinL^{p) 



< 



* ds 



I n^, ,l[\\m+knflis)J2{mL^iL^nL^) + \\nhis)j2is)]h.Lr 

Iq (1 + r — Sj V q \ X • n; I q l 

When £ > 3/6 — 1 (since fc < 1), using interpolation the above is bounded by 

ds 

Next Lemma [Ql allows us to bound the above by 

^ I (l+f-g)fc ll™^/iWllLr(i^nL~) \M2{mL^iLlnL^) ■ 
From Proposition 13.51 and p.l2p we see that the last line is 

/■* ds 

< ||n7fcu;£/i||i»=^(i2nL~)||rofcW£/2||L- (L2nL~) (1 + f - s)fc(l + s)2fc 

< (1 + ty''\\wkWefi\\L^^(^LlnL^)\\'^kWef2\\L^^^(LlnL^)- 
We have shown 

||n7fcWf7V[/i,/2]||L~ (L2nL~) < \WkWefi\\L^^(LlnL-^)\WkWef2\\L-^^{LlnL^)- 

To handle the linear semigroup, U{t), wc again use Theorem 13.11 and Theorem 14.11 
to obtain 

< (l|w£+fe/o|lL~(L2ni^) + WfahlL- + W^kWtfWl^^^LlnL^)) ■ 

We conclude that M[-] maps M^^ into itself for < i? chosen sufficiently small and 
e.g. ||w£+fe/o||2^cx,(2,2nL°=) + ll/ollL2 2,r < 2^77' obtain a contraction, we consider 
the difi^erence 

A/[/i] - A/[/2] - N[f, - /2, /l] + iV[/2, /l - /2]. 
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Then as in the previous estimates we have 

\\rv,w,{M[h]~M[m\L^^^LlnL^) 

< Clk (ll^fe«^f/l||L~ (L2nL~) + \ \^kWt.f2\\L^JLlnL^)^ 

X \\mkWg{fi - f2)\\L^jLlnL^)- 

With these estimates, the existence and uniqueness of sohitions to (|2.5p follows 
from the contraction mapping principle on Mj^^ when i? > is suitably small. 
Step 2. Continuity. We perform the estimates from Step 1 on the space 

= C'iiO, oo) X X Rl) n M^, R>0. 

As in Step 1, we have a uniform in time contraction mapping on M^f for suitable 

R. Furthermore M[f] is continuous if / G ^^k f ^'^'^ /o continuous. Since the 
convergence is uniform, the limit will be continuous globally in time. This argument 
is standard and we refer for instance to [21][24l|3T] for full details. 

Step 3. Positivity. We use the standard alternative approximating formula 

with the same initial conditions -P'"^^|(_q = = ^ + v^foi for n > 1 and for 

instance =^ J + VJ/o- Here we have used the standard decomposition of the 
collision operator Q = — Q_ into gain and loss terms with 

Q_(F"+i,^^") = / dq f duv^ a{g,0) F"+\p)F'\q) = 

and = Q_(1,F"). If we consider F''+^{t,x,p) = J + Vjp+\t,x,p), then 

related to Step 1 we may show that it;£(p)/"+^(t, x,p) is convergent in L^((L^nL!^) 
on a local time interval [0,T] where T will generally depend upon the size of the 
initial data. In particular f"~^^{t,x,p) = ^ satisfies the equation 

{dt+p-v^ + v{p)) .r+i = Kin + r+(.r, n - r^{.r+\n. 

Here again we use the standard definition r+ — r_ = F where 

T_{r + \r)= [ dq f dcov,aig,e)^/ji^) r+\p)r{q). 
We rewrite this equation using the solution formula to the system p.2p as 

This solution formula G{t) is defined just below p.2p . Furthermore 

cir^\,r) = f ds G{t - s)K{r) 

Jo 

+ [ ds Git - s)F+(/", /") - Git - s)F_(/"+\ /"). 

"'0 

For given T > and i? > we consider the space M^^([0, T]) defined by 
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We can do something precisely similar to Step 1 to prove the existence of /"^^ E 
M^^([0,T]). Instead of using Theorem [O and Theorem Ol to deal with the 
operator [/, we use in particular Lemma 13.21 in this paper and Lemma 4.2 in |46| 
to deal with the operator G. Then we have estimates similar to those in Step 
1 for our mapping in the space M^^([0,r]). Now given /" e A/^(,([0,T]) and 
|jwf+fe/o||^oo(-^2p,2.°°) + ll/o|lL2^r < I with R > chosen sufficiently small, as in 
Step 1, we can prove the existence of /"+^ e ^fe^f ([0, T]). 

With the estimates established in this paper, it is now not hard to show that 

Here T > is sufficiently small, and the constant C > can be chosen independent 
of any small T. Therefore there exists a T* > such that w^f^ — > w^f uniformly 
in L^{L^ n L^) on [0,r*]. This will be sufficient to prove the positivity globally 
in time. 

Indeed if i^" > 0, then so is F^) > 0. With the representation formula 

pt 

+ / ds e-J'=''^"^(^"'("'"-P(*-")'P) Q+(F",i^")(s,a;-p(t-s),p). 
Jo 

Induction shows F"^^ {t, x,p) > for all 7i > if Fq > 0, which implies in the limit 
n — > (X) that F{t,x,p) = J + ^/j f{t,x,p) > 0. Using our uniqueness, this is the 
same F as the one from Step 1 on the time interval [0, T*]. We extend this positivity 
for all time intervals [0,T*] + T*k for any fc > 1 by repeating this procedure and 
using the global uniform bound in from Step 1. □ 
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